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Abstract — This article discusses the steady analysis of an 
M/M+G/2 queue. All arriving customers are served either by 
server-1 according to an exponential service time density 
f‘i(x)=pe'' ,s with mean rate p or by server-2 with a general 
service time distribution B(t) or density function b(t)dt=dB(t) 
with mean rate p 2 =l/|i or mean service time is p. Sequel to some 
objections raised on the use of the classical 'First Come First 
Served (FCFS)' queue discipline when the two heterogeneous 
servers operate as parallel service providers, alternative queue 
disciplines in a serial configuration of servers are considered in 
this work; the objective is that if, in a single-channel queue in 
equilibrium, the service rate suddenly increases and exceeds the 
present service capacity, install a new channel to work serially 
with the first channel as suggested by Krishnamoorthy (1968). 
Using the embedded method subject to different service time 
distributions we present an exact analysis for finding the 
‘Probability generating Function (PGF)’of steady state number 
of customers in the system and most importantly, the actual 
waiting time expectation of customers in the system. This work 
shows that one can obtain all stationery probabilities and other 
vital measures for this queue under certain additional and 
simple but realistic assumptions. (Abstract) 

Index Terms — Poisson arrival, service time distribution; 
PGF of queue length distribution, waitime distribution, mean 
queue length and mean waiting time 


I. Introduction 

We analyse a class of M/G/2 queueing models to study the 
needs for designing serially connected system with two 
heterogeneous servers/machines which does not violate the 
FCFS principle in place of paralleled heterogeneous servers 
which will violate the classical FCFS queue discipline. 

For instance, if there are two parallel clerks in a reservation 
counter who provide service with varying speeds then 
customers might prefer to choose the fastest among the free 
servers, other hand if the customer chooses the slowest server 
among the free servers then the one who could come 
subsequently may clear out earlier by obtaining service from 
the server providing service with faster rates what is exactly 
called a violation of the FCFS discipline. 

Hence there is a need for the designing of alternative ways to 
the parallel service providers that would reduce the impacts of 
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violating the FCFS so that the resulting waiting times of 
customers are identical. 

Customers arrive according to a Poisson process at a rate X for 
service on these servers with an expectation of spending least 
waiting time. 

For customers serviced by server- 1, their service time Ti 
follows exponential distribution with rate p i.e. 
Fi(t)=P(T !<t))= 1 -e _fit and the probability density 
function(PDF) is fi(t)=dFi(t) and hence the ‘Laplace 
Transform(LT)’ of fi(t) is fi (s)= 

e ' (l )dt = — — — ; Re s>0 . Let the service times T 2 of 
s + JU 

customers serviced by server-2 follow a common 
distribution B(t) with PDF b(t) and mean p and the LT of b(t) 

OO 

is b*(s) = j”e "b(t)clt . Both of these service time 

o 

distributions are assumed to be mutually independent and 
each is independent of inter-arrival time distribution also. 

A simple way of connecting the servers in series subject to 
servicing of customers according to the FCFS queue 
discipline is proposed below with three types of service time 
distributions. For a single server queue of M/G/l type under 
equilibrium conditions, expecting greater demand for service 
in the long run it is decided to increase the service capacity of 
the primary single server- 1 by installing an additional 
server-2 as detailed below: 

• If an arriving customer enters into the idle system, his 

service is immediately initiated by server-1. This 
customer is then served by the server- 1 at an exponential 
rate p if no other customer arrives during the on-going 
service period; 

• Otherwise i.e. if at least one more customer arrives before 

the on-going service is completed then the same customer 
is served jointly by both servers according to the service 
time distribution F min (t)=P(T min <t), where T min = Min(Ti, 
T 2 ) and the PDF of T rllln is f mm (t). 

Thus if system size N(t) i.e. the number of customers present 
in the system at time ‘t’ is greater than or equal to 2, then 
server-2 joins server- 1 to server the customer according to the 
service time distribution F^nCO; otherwise the server- 1 is 
alone available at the service counter. The LT of F^lt) is 
00 

( 0 * 

0 

=f;( S )+b\s+ J u){i-f 1 \s)} 

It is attempted here to study the two-server (heterogeneous) 
M/M+G/2 queues in the light of the above queue discipline 
with different service time distributions F min (t), F max (t) and 
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Fmin+max(t:) while the servers are connected in a series system. 
We are motivated by the numerous physical applications and 
contributions since the above M/M+G/2 model can 
conceptualize well in the analysis of queues found in modern 
telecommunications, computer business centres, banking 
systems and similar service operations systems where human 
nature is evident. Over the years, a lot have been written on 
homogeneous service systems as in Hoksad [3], The reader is 
referred to and [4] through [8] and many others to refresh on 
the heterogeneous service systems. A model of the general 
service type is studied by Boxma, Deng and Zwart [2], 
unfortunately, due to complex structuring, formations and 
assumptions; it could not estimate certain areas in the general 
case. 

The probability generating function (PGF) of the number of 
customers present in the system, LT of the waiting time 
distribution and their mean values have been obtained in 
section II. A numerical illustration is also provided to support 
the results on mean waiting times. Subsection C highlights 
the various special features of the proposed methodology and 
its future scope. 


II. STEADY STATE CHARACTERISTICS OF M/M+G/2 
QUEUES WITH SERVICE TIME DISTRIBUTION 

The number N(t) of customers present in the system at time t, 
does not now constitute a Markov process, see [1], We 
consider here 'the embedded time points' generated at the 
departure instants of customers just after a service completion 
either by server- 1 or by server-2. We consider the Markov 
Chain of system states at these embedded points where the 
state of the system is represented by the number, N k =N(t k ), of 
customers left behind in the queue by the k th departing 
customer at the departure epoch ‘t k ’. The discrete time 
process {N k } constitutes a Markov chain on the state space 
S ={0,1,... oo }. 


A. PGF of {Nk} 

Let qj be the steady state probability of finding 'j' customers in 
the system as observed by a departing customer and the 
z-transform of the probability distribution of {q, ; j=0, 1, ,... 


oo 


be V(z) = ^ j q j Z j . 

j = o 

Let a, denote the probability of 'j' arrivals in a service time 
PDF f m in(x) and 5j denote the probability of 'j' arrivals in the 
exponential service time PDF fi(x). Since the arrivals come 
from a Poisson process at rate L, we get, for j=2, 3, ... oothat 


a > = 


and 




./'! 


■f m Jx) dx ( 1 ) 


—Ax 


(Ax) J 


J 


i 1 


fi(x) dx 


( 2 ) 


Let the respective z-transforms of the probability distributions 
{dj} and {5j } be A min (z) and A^z): 

00 

A m „(^) = Z a i ZV =/min*0*-^)(3) 

7=0 

and 


A (z) = Yj 5 J A = fi(A-Az) (4) 

7=0 

Focusing on the embedded points under equilibrium 
conditions, let the unit step conditional transition probability 
of the system going from state 'i' of the tk-l f embedded point 
to state 'j' in the k th embedded point be qjj=P(N k = j /N k _i=i) 
for i, j e S. These transition probabilities will form the unit 
step transition probability matrix Q=f qy) as below: 


Q 


(8 

u 0 

A 

A 

A 



A 

A 

A 

A 



0 

a 0 

a l 

a 2 

a 3 


0 

0 

a 0 

a x 

a 2 

a 3 


0 

0 

a 0 

«[ 

OC 2 OL 3 -y 


(5) 


Thus equilibrium state probabilities at the departure instants 
are given by qj = lim q j " where q represents the n-step 

n — >oo J 

probability of moving from state 'i to j'. Let q=(q 0 ,q] ,q 2 ,. . .) be 

arow vector and let e=(l,l, , 1) be a column vector of unit 

elements of infinite order. Assume that A nlin ’(l) < l, then the 
stationary distribution of the state transition matrix Q, exits 
and is given by the unique solution of the following system of 
equations: 

q Q = q and q e = 1 (6) 

Multiplying the j lh (j=0, 1,2, . . .) equation of qQ=q , of (6) by 
z J and summing all the left-hand sides and the right-hand 
sides from j=0 to j= co , we get the PGF(Probability 
generating Function) V Ill]n (z) of the queue length distribution 
{qj} of the sequence{N k }. 

V ■ (z) = 

min \**s 

glZ[Anin(z) ~ A(Z)1 + gptAmnU) ~ A(Z)Z] ^ 

An „(Z)“Z 

1 X 

Let JU 2 and p — — be used in the sequel. It is remarked 

- P p 

that closed form expression to b (s+p) cannot be obtained 
unless particular cases like LT of exponential, Erlangian, 
Phase (PH) type or hyper exponential distributions. Assume 
that the service time distribution of the server-2 is one of these 
distributions such that b (p) is finite as a function of p and p 2 
for our discussions to follow. Since q t = pq 0 , A lmn '( 1 )< 1 . 
A['(l)=p and V(l)=l, we derive from (7) that 


qo : 


l-^And) 


( 8 ) 


lHp-A'^mx+p) 

Thus v min (z)= 

(Z ~ 1) A (z)] + (1 + pz)[ A (z) - Anm (z)] 


Z-Ann(z) 


l-^'nnnd) 


(9) 


l+Go-AvoMi+P) 

The mean number E(N) of customers present in the system at 
a random point or at a departure epoch of time is E(N)= 

'{p-A^'mp+A^'m+p)] 


+ p+ 


! + {/>" Am '(1)}(1 + P) 

■Anin '(l) 2 


1-A.n'd) 


( 10 ) 
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Denote the LT of the waiting time W distribution 
W(t)=P(W\(< t) by w(s) for 0<s<l/A... Then replacing z by 
(1-s/A.) in (9) and using the fact w (s)=V(l-s/A.), we derive 
w (s) as 


l- A' 


w (s) = 


.( 1 ) 


i hp-a' . (xm+p) 

min 


(11) 


{ 1 + pil - s / X) } { Af . * (s) - Af .* ( s) } + sf .* (s) 


s-A{\-f . (^)} 

nun 


The mean waiting time W of a customer in the system 
obtained from (1 1) is found to satisfy the well-known Little's 


formula X W =E(N). 


B. Design of Shortest Delay Time Environment for 
M/M+G/2 Queues 


The forgoing embedded methodology of analysing the 
M/M+G/2 queues can also be extended to design a shortest 
processing environment as discussed below. 

Let F max (t)=P(T max <t), where T max = Max(T,, T 2 ) and the PDF 
of T max is f max (t). The LT of f max (t) is 
00 

/«■*/_ (0* 

0 


= b\s)-b\s + fi){l-f l \s)} (12) 

Considered is the joint service time PDF ‘f(t)’ of the two 
servers 

f(t) = 71 1 f mi n(t) + Jt 2 f max (t) (13) 

in place of f min (t) where rcj and ;t 2 are the probability values 
of choosing the service time PDFs f min (t) and f max (t) 
respectively. 

Let the corresponding z-transform of number of arrivals 
during the joint service PDF f(t) be 

oo 

A{z) = ^a j z i = x l f rt A{A-Az) 

j=o 

+ n 2 f m J'(A-Az) 

A /l 

Let p 2 = — and p— — . Assume that A’(l)< 1 of(14), i.e. 

th P 


A\\) = x l p{\-b*( 1 u)} 

, > (15) 

+ x 2 {p 2 +pb*(p)}<l 

Now on substituting A(z) of (14) in place of A lnill (z) of (9), 
one can verify that (V(z) replaces V^lz) ) 


V(z) = 


(z-l)4(z) + (l + /Pz)[A(z)-A(z)j 
z-A(z) 


l-A’(l) 

l + (/?-A'(l))(l + /?) 


(16) 


As before, to get the closed form expression to the 
w*(s)=V(l -s/a) i.e. LT of the waiting time distribution, it is 
obvious that the factor z must be replaced by ( 1 -sA.) on both 
sides of (16). Thus both the queue length and the waiting time 
distributions of the M/M+G/2 queues are completely 
determined even if the joint service time PDF of the two 
servers is f(t) of (13). 


Lemma : If p 2 and n x < 7T-, then A’(1)<1 of (15) 


implies that p + p^ < — . 

n x 

Let W| and w 2 be the mean waiting times corresponding to p < 
(j, 2 and P >!-i ? values respectively. Thus w 1 =E(N)A when p < 
p 2 and w 2 =E(N)/k when p > p 2 and 

A\lf 

P + — + 


E(N)= < 


l-A'(l) 

{p-A\\)}[p + A'm + p)\ 

i+{p-Axmi+p) 


(17) 


From the design criterion of the above model, it is expected 
that if the service rate p of server- 1 is larger than that p 2 of 
server-2 then w 2 < Wj with increasing values of Jp. For an 
illustration, w t is calculated fixing /,~4.8, p=4.5 and p 2 =8 
(exponential case) and w 2 is computed fixing k=4.8, p=8 and 
p 2 =4.5 while Jti varies between 0.4 and 1. Fig 1 below shows 
the graph of these variations on 1 OOOw and 1000w 2 which 
also ensures that w 2 <W) uniformly. 


2.5 



Fig 1: Graph showing tii values versus mean waiting times w, 
and w 2 . 


This is a kind of unequal service rates queueing application 
where the event of shortest mean delay happens only if the 
service rate p of server- 1 is larger than that p 2 of server-2. 


III. Conclusion 

As the analysis of the above M/M+G/2 queues seems to be 
simpler and it is an easy procedure to implement if any real 
life application warrants. We conclude that the proposed 
method so far discussed is a good alternative for use instead of 
fitting of M/G/2 queues. Most importantly there are few new 
results produced from above work that could be applied on 
the two-serially connected machines in the field of 
engineering problems rather than applied science areas. To 
extract more information like arrival epoch probabilities one 
can employ Markov Renewal theory as in Senthamarikannan 
and Sivasamy [9] and [10]. 
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